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Abstract 

Remarkable progress inspired by twistors has lead to very simple analytic expressions 
and to new recursive relations for multi-parton color-ordered amplitudes. We show 
how such relations can be extended to include color and present the corresponding 
color-dressed formulation for the Berends-Giele, BCF and a new kind of CSW recursive 
relations. A detailed comparison of the numerical efficiency of the different approaches 
to the calculation of multi-parton cross sections is performed. 
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1 Introduction 



The quest for physics beyond the standard model will soon enter an exciting phase as the LHC 
starts colliding protons at 14 TeV. At such an energy we will be probing the mechanism which 
breaks the electroweak symmetry and hopefully understand what stabilizes the associated 
scale. Whether this mechanism can be described perturbatively by a standard model Higgs 
or can be related to new symmetries (such as supersymmetry) or even to the discovery of an 
extended space-time structure, is the subject of an intense theoretical activity. The possibility, 
however, of answering such far reaching questions at the LHC will eventually depend on our 
ability to discriminate signals of new physics from large standard model backgrounds. With 
so much energy available in the center of mass, the cross section for multi-jet QCD events 
even associated with heavy objects such as weak bosons or top quarks will be in most cases 
very large if not overwhelming. 

As the need for better predictions for QCD processes with many particles in the final state 
has become clear, a substantial activity on developing the techniques and the tools to perform 
such calculations has spurred. Several codes have by now become available that can compute, 
numerically, tree-level cross sections and generate events with many particles in the final state 
in an automatic way [1-4]. Progress has been significant also in improving the accuracy of the 
predictions by calculating next-to- and next-to-next-to-leading QCD corrections for processes 
with up to three and just one particle in the final state respectively [5]. 

In the midst of this effort, unexpected theoretical progress has come from the so-called 
twistor-inspired methods [6], which have provided new techniques to compute analytic results 
for gauge amplitudes both at tree and one- loop level [7]. These new methods go back to a 
correspondence between a weakly coupled M = 4 Super Yang-Mills theory and a certain type 
of string theory. The key point in this correspondence is that all tree-level (color-ordered) 
amplitudes are related to algebraic curves in twistor space. In Ref. [8] it was shown that this 
leads to the so-called CSW rules, which state that all tree-level (color-ordered) amplitudes 
can be constructed using a small class of very special amplitudes, the maximally helicity 
violating (MHV) amplitudes. Another important result has been the formulation of a new 
kind of recursive relations in addition to the well-known Berends-Giele recursion [9], the so- 
called BCF relations [10,11]. They state that any tree-level (color-ordered) amplitude can be 
constructed from products of two on-shell amplitudes of fewer particles, multiplied by a simple 
scalar propagator. These new calculational tools have allowed to derive expressions for some 
multi-parton amplitudes [8,12,13] and certain classes of splitting functions [14,15], which have 
simple and compact analytic forms. So far, many extensions of the twistor-inspired methods 
have been presented, in particular generalizations to include scalars [16], fermions [17,18] and 
photons [19]. 

The purpose of this work is to present a new extension of the twistor-inspired methods. 
We show that it is possible to reformulate the twistor-inspired recursive relations for the 
color-ordered amplitudes in terms of the full amplitudes. The motivations are twofold. The 
first is mostly theoretical and stems from the observation of an interesting similarity between 
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a color decomposition based on the adjoint representation and the BCF recursive relations 
which suggested the existence of a formulation that embodies both. The second is more 
pragmatic and aims at establishing whether the new twistor-inspired recursive relations are 
an improvement also at the numerical level. In order to make a consistent comparison with 
the most efficient algorithms available [1,20], a recursive formulation which includes color 
is necessary. In the standard approach where the color-ordered amplitudes are calculated 
analytically (or numerically), one has to sum over the permutations of the color orderings 
to obtain the full amplitude. This algorithm has an intrinsic factorial growth and cannot 
compete with the available numerical methods which only grow exponentially. In this work 
we derive a general method to reinstate color into the recursive relations for color-ordered 
amplitudes and apply it to the BCF and to a modified version of the CSW relations. 

The paper is organized as follows. In Section[2]we review the notion of color decompositions, 
highlighting the different features of the various color bases available. In Section El we discuss 
the Berends-Giele recursive relations and we present a simple derivation of their color-dressed 
counterpart which serves as an illustration of the method that will be applied later. In 
Section |U we prove our first main result, i.e., the color-dressed version of the BCF relations, 
Eq. (|4.12jl and discuss its most important features. In Section El we reformulate the CSW 
relations in terms of simple new three-point effective vertices and derive their color-dressed 
version. Section El contains the numerical results on the evaluation of multi-gluon amplitudes 
obtained using the different color-dressed recursive relations, focusing on the comparison with 
known techniques. Finally we draw our conclusions. 

2 Color Decompositions 

In this section we briefly review the notion of color decomposition of tree-level QCD ampli- 
tudes and the available results. Emphasis is given to those aspects that will play an important 
role in the following. 

The basic idea of a color decomposition is to factorize the information on the gauge struc- 
ture from the kinematics. As an example, consider the amplitude for n gluons of colors 
Oi, 02> • • • ? a n ( a i — 1) • • • ) N 2 — 1)- O ne can easily prove that at tree level, such an amplitude 
can be decomposed as [21] 



where T a are the fundamental-representation matrices of SU(N), and the sum is over all 
(n — 1)! permutations of (2, . . . , n). Each trace corresponds to a particular color structure. 
The factor associated with each color structure, A, is called a color-ordered amplitude. 1 It 
depends on the four-momenta Pi and polarization vectors of the n gluons, represented 
simply by i in its argument. The color-ordered amplitudes are far simpler to calculate than 

^^Also referred to as a dual amplitude or partial amplitude. 
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the full amplitude A due to the smaller number of Feynman diagrams contributing to them. 
They have several remarkable properties. Among them, a special role is played by the so- 
called Kleiss-Kuijf relations [22]. These are linear relations amongst the amplitudes directly 
inherited from the gauge structure, i.e., from color, which in the case of n-gluon amplitudes 
reduce the number of linearly-independent amplitudes to (n — 2)!. It is then clear that the 
number of terms in Eq. (|2.1j) is not minimal. This decomposition has no special feature except 
that it was inspired by string theories (color factors are the Chan-Paton factors of the open 
strings). It is universally used to illustrate the idea of color decomposition and to define the 
color-ordered amplitudes A in terms of the full amplitude A. It can be shown, however, that 
this definition does not depend on the color basis. 

Recently, another decomposition has been introduced, which is based on color flows [23,24]. 
This decomposition arises when treating the SU(N) gluon field as an N x N matrix (A^j 
(i,j = 1,...,N), rather than as a one-index field A ^ (a = 1,...,N 2 — 1). The n-gluon 
amplitude may then be decomposed as 

A(l,...,n)= Yl C 2 C 2 3 --^t^(l>^...,^), (2.2) 

where the sum is over all (n — 1)! permutations of (2, . . . ,n). The partial amplitudes that 
appear in this decomposition are the same as in the decomposition in the fundamental rep- 
resentation. The color-flow decomposition has several nice features, the most attractive one 
being its simplicity. In fact, it is a very natural way to decompose a QCD amplitude, and as 
the name suggests, it is based on the flow of color, so the decomposition has a simple physical 
interpretation. It does not involve any color matrix and the color factors in front of each 
amplitude are either zero or one. This basis makes the numerical calculation of color factors 
in the evaluation of full amplitudes fast, even though the number of terms in the sum over the 
(n — 1)! permutations is not minimal. Finally, such a decomposition exists for all tree-level 
parton amplitudes including any number of quark pairs and gluons. 

A third decomposition of the multi-gluon amplitude is available, which is based on the 
adjoint representation of SU(N) rather than the fundamental representation [25,26]. The 
n-gluon amplitude in this decomposition may be written as 

A(l,...,n)= (F a °>F a °>--.F a -^) aian A(l,a 2 ,...,a n ^,n) , (2.3) 

where (F a ) b c = —if abc are the adjoint-representation matrices of SU(N) (f abc are the structure 
constants), and the sum is over all (n — 2)! permutations of (2, ...,n — 1). The indices 
corresponding to the first and the last gluon are taken as "references" and are not included in 
the permutations. The partial amplitudes that appear in this decomposition are the same as 
in the other decomposition, but only the (n — 2)! linearly- independent amplitudes are needed. 
In this respect this formulation is "minimal" as there is no redundancy and the Kleiss-Kuijf 
relations are embodied in the color factors. As we will elaborate upon in the following, there 
exists a remarkable formal similarity with the BCF recursive relations, where two gluons are 
also taken as a reference to build up the full amplitude. 
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3 Color-dressed Berends-Giele relations 



In Ref. [9], Berends and Giele introduced the color-ordered n-point gluon off-shell current 
J M , which can be defined as the sum of all color-ordered Feynman diagrams with n external 
on-shell legs and a single off-shell leg with polarization \i. The color-ordered off-shell currents 
can be constructed using the Berends-Giele recursive relations 



{n-l 
(Pi,k, Pk + i,n) J*(l, ■ • • > k)J p (k + 1, . . . ,n) (3.1) 
k=l 

n— 2 n— 1 

+ E E Vj a ""Ul,...,j)J P (j + l,...,k)Uk+l,...,n) 

3=1 k=j+l 

where 

Pi,j =Pi+ Pi+i + • • • + Pj-i + Pj, (3.2) 

and V^ vp {Px t k, Pk+i,n) an d V± vpa are the color-ordered three and four-gluon vertices defined 
in Ref. [27]. It is easy to see that the four-gluon vertex appearing in these relations introduces 
a larger number of possible combinations of subcurrents than the three-gluon vertex. It is 
possible to simplify the recursion by decomposing all four-gluon vertices into three-vertices 
including a tensor particle (Fig. Using this decomposition, the Berends-Giele recursive 
relations can be rewritten such that only three-point vertices are present 

. n-l ( 

J"(l, 2, . . . ,n) = - J- \ v r ( p i,k, Pk + i,n) Ml, ■ ■ ■ , k) Mk + 1, • ■ • , n) (3.3) 
i'™ k=i { 

+ V^Ml, k)J a p(k + 1, . . . , n) + V^J aP {l, . . ., k)J a (k + 1, . . . 
where J a p is a tensor off-shell current, and Vp ua/3 is the tensor-gluon vertex, defined as 

vr a = f isTsT - rg vp ) ■ (3.4) 

As there exists no one-point tensor off-shell current, all such currents appearing in Eq. ()3.3j) 
are defined as zero. The tensor off-shell currents can be easily constructed recursively from 
gluon off- shell currents 

n-l 

V(l,2,...,n) =iD^ V° paP J2M^--^k)Mk + l,---,n), (3.5) 

k=l 

where iD^ap is the color-ordered tensor "propagator" , defined as 

iD lu , pa = - t -{gPtf°- ! nr). (3.6) 
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Fig. 1: Diagrammatic representation of the decomposition of the color-ordered four- 
gluon vertex. 

We now present a systematic method to dress color-ordered recursive relations with color 
in order to obtain recursive relations for the color-dressed off-shell currents. In the color-flow 
decomposition, a color-dressed gluon off-shell current can be written as 

Jtfl, 2, . . . ,n) = 5 L S tl ■ ■ ■ flf^K °2, ■ ■ ■ , ( 3 - 7 ) 

where (/, J) is the color of the off-shell leg. A color-dressed tensor off-shell current can be 
obtained similarly. We will explain the color dressing of the Berends-Giele recursive relations, 
Eq. (J3.3J) . dealing with the pure gluon vertices and the tensor-gluon vertices separately. After 
inserting Eq. ()3.3j) . into the color-flow decomposition, Eq. fl3.7|) . the three-gluon vertex part 
reads 

n-l 

^EE 6 L 6 tl ■ ■ ■ ^ nV * VP {. P ^ P ^n) -, cr k )J P (a k+1 , ...,a n ), (3.8) 

1 ' n a£Sn k=i 

where 

P(Ti ,cr k = Pa\ ~\~ Pa 2 + • • • + Pa k j 

P(rk+l,<*n ~ P°k + 1 P<?k+2 + • • • + Pa n 

The color factor appearing in Eq. ()3.8|) can be written as (Fig. |2J 

5{jc 2 ...^r = 5 J G 5f sgs§6§ C fc+1 ■■■sir, (3.10) 

where 

- 8 G 8f is the color structure of the propagator appearing in the Berends-Giele recursive 
relations. 

- 8f . . . 8 J ^ k is the color structure of the subcurrent J v {p\, ■ ■ ■ , <?k), where the off-shell leg 
v has color (K,L). 

- 8^ i . . . 8 3 ^ is the color structure of the subcurrent J p ((Tk+i, ■ ■ ■ , cr n ), where the off-shell 
leg p has color (M, N) . 

5 



(3.9) 




Fig. 2: Diagrammatic representation of the decomposition 1)3. lOj) of the color factor 
of the three-gluon vertex part. 

- 5 ^ 5$ tiff is part of the color structure of a three-gluon vertex to which the off-shell legs 
fi, v, p with colors (G, H), (K, L), (M, N) are attached. 

We now define an ordered partition of a set E into two independent parts as a pair (7Ti, %%) 
of subsets of E such that 7Ti © 7r 2 = E, which means (7Ti, 7r 2 ) 7^ (7r 2 , 7Ti). Furthermore, we call 
(unordered) partition of a set E into two independent parts a set {7ri,7T2} of subsets of E 
such that 7Ti © 7r 2 = .E and {7Ti, 7r 2 } = {7r 2 , 7Ti} . These definitions can be easily extended to 
partitions of a set E into n > 2 independent parts, for both the ordered and the unordered 
case. 

In the case encountered here E = {1,2, ... ,n}. We will denote the set of all ordered 
partitions of E into two independent parts by OP(n, 2) and the set of all (unordered) partitions 
of E into two independent parts by P(n, 2). Using these definitions, the sum over permutations 
appearing in Eq. (J3.8)) can be decomposed as follows: For a given value of k, 

- Choose an ordered partition n = (7Ti, 7r 2 ) in OP{n, 2) such that #71"! = k, where j^T\\ is 
the number of elements in the set ~R\. 

- Fix the first k elements of the permutation to be in the subset tti. 

- Sum over all permutations of the first k elements and over all permutations of the last 
n — k elements. 

- Sum over all possible choices for the ordered partition % = (jn, 7r 2 ). 
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This is equivalent to the replacement 



ra-1 



E E - E EE- (3-m 

fc=l crG5„ 7TGOP(n,2) cre5 fe cr'eS n _ fe 



The three-gluon vertex part now reads 

E E E ^ss§vr(p w ,p Kt+lK ,] 

1 ' n 7rgOP(n,2) <TGS fe <r'G5 n _ fc 

. . 4 %fe J„(^ , . . . , a n u) 6* . . . 5^" J P « k+ i, ■■■,<«), (3.12) 
where 7i"i = {7T 1 , 7r 2 . . . , 7r fc } and 7T 2 = {7r fe+1 , 7r fc+2 , . . . , 7r™}. 

Clearly, k and P CT ' )(T ' „ only depend on the choice of the ordered partition n = (ttx, 7r 2 ), 

but not on the order of the elements in tt\ and i\2- We therefore define 



(3.13) 



-Ptt 2 = Ptt^ 1 + PV+ 2 + • • • + Pn n - 

It is now possible to identify several subcurrents in this expression, namely 

jf L M = i.-^k-,^ (3.i4) 

^fo) = E C (3-15) 
so that the three-gluon vertex part reads 

Wjj- E j v kl mj?*m- (3.i6) 

1,n 7reOP(n,2) 

In Ref. [24] it was shown that the (color-dressed) three-gluon vertex can be expressed in the 
color-flow decomposition as 2 

vr (p n ,p n ) = fis&g vr (p^p^) + 4sM vr (p^p.j . (3.17) 

So, finally the three-gluon vertex part can be written as 

t&fj^- E v$ vp {P*i>P*>)tf L MJi i *M. (3-18) 

1 ' n 7rSP(n,2) 

We now address the tensor-gluon part in the Berends-Giele recursive relations, Eq. (J3.3|) . 
From Fig. El it can be seen that for the s-channel appearing in the decomposition of the 
four-gluon vertex, each tensor-gluon vertex has the same color-flow structure as a three-gluon 
vertex. The t-channel contribution is similar. Therefore the tensor-gluon vertex can be written 
in the color-flow decomposition as 

2 For brevity, the color indices of the vertex are not written explicitly. 
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Fig. 3: Color-flow structure of the s-channel contribution to the for gluon vertex. 



V. 



x J S3i XJ2 i/VvaP , r J rj2 srji i/^mf3 
h°i2 I V T ^ %2°i\ l V T ) 



(3.19) 



where (J, J) is the color of the tensor particle. The color dressing of the tensor-gluon part in 
Eq. ()3.3|) is hence exactly the same as for the pure gluon part, leading to 

E { v^j^mj^m + vr'Jaf^jr^) } . (3.20) 

7r€OP(n,2) 

As the sum runs over all elements in OP(n, 2), we can exchange 7Ti and 7r2 as well as the color 
indices (K, L) and (M, N) in the last term. Using Eq. 1)3.19)1 the tensor part now becomes 



E vt^j^mj^w 

7rgOP(n,2) 



(3.21) 



Hence the color-dressed recursive relations with all four-gluon vertices replaced by tensor 
particles read 



J, 



i J i 



D^tj (Pl,n) 



+ 



7reP(n,2) 



E vt^j^mj^m 

7rSOP(n,2) 



(3.22) 



To complete the color dressing of the Berends-Giele recursive relations, we have to apply 
the color-dressing method introduced above to the recursive relations for the off-shell tensor- 
currents, Eq. (J3.5j) . Both the color-dressed vertex and the color-dressed off-shell tensor-current 
have the same form as in the pure gluon case and the recursive relations for the tensor particle, 
Eq. (|3.5|) . have the same structure as for the three-gluon vertex part in the previous section. 
Therefore, one can immediately write down the color-dressed recursive relations for the off- 
shell tensor-current 



n) 



(3.23) 



7rGP(n,2) 

The two recursive relations, Eq. ()3.22|) and Eq. ()3.23|) . can be solved simultaneously to con- 
struct color-dressed gluon off-shell currents for arbitrary n. The full color-dressed scattering 



S 



amplitude is then recovered be putting the off-shell leg on-shell. This result is equivalent to 
the Dyson-Schwinger algorithm presented in Ref. [4]. It should be noticed that the color- 
dressed recursive relations have the same form as the color-ordered Berends-Giele recursive 
relations, Eq. (jH.HJ) and Eq. (|3.5jl . The only difference between the color-ordered and the 
color-dressed case is that in the latter we sum over unordered objects and no permutations 
need to be taken into account. This is a general feature which will turn out to be common to 
all color-dressed recursive relations. 



4 Color-dressed BCF relations 

In this section, we apply the same method employed to construct the color-dressed Berends- 
Giele recursive relations to the BCF recursive relations, presented in Ref. [10,11]. Assuming 
that gluons 1 and n have opposite helicities, the BCF recursive relations read 

n-2 ^ 

A n (l, 2, ...,n) = Yl A ^ 2 > • • • > fc ' ~ P i~k) JxT A n-k+i (AV k+l,...,nj , (4.1) 

k=2 1 >* 



where a sum over the helicities h of the intermediate gluon is implicit, and 

P 2 

Pi,k = Pi,k+ , ip'^A Ai, (4.2) 
P 2 

Pi = Pi + / i p' fc hi A " A i' ( 4 - 3 ) 
P 2 

Pn = Pn~ / p'^i Ul, (4.4) 

with Aj and Aj being the spinor components of pi = Xi Aj. 

As in Eq. (14. lj) we have to choose two reference gluons, 1 and n, the color-flow decomposition 
and the color decomposition in the fundamental representation are not well suited to dress 
the BCF recursive relations with color, because they allow us to fix only one of the two 
reference gluons. The most natural color decomposition which fixes both reference gluons is 
the color decomposition in the adjoint representation, Eq. (J2.3|) . Inserting the color-ordered 
BCF relations, Eq. IfP]) . into Eq. (J2~3j) . one finds 



n-2 



An(l,...,n) = J2 E (^•••^-- 1 ) aian ^ + i(i^2,...,^,-P 1 



k=2 aGS n -2 



— \— A„_ fc+ i (P\^ k , (Tk+i, • • • , cr n -i, nj , (4.5) 

where 

Pl,a k = Pi + Pa 2 + • • • + Pa k ■ (4.6) 
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For a given value of k, the sum over permutations appearing in Eq. (|4.5|) can be decomposed 
in a similar way as for the three-gluon vertex part in the Berends-Giele recursive relations. 
The procedure is as follows 

- Choose an ordered partition it = (tti, 7^) of {2, 3, . . . , n — 2, n — 1} such that #7Ti = A;— 1. 

- Fix the first k — 1 elements of the permutation to be in the subset 7!"i. 

- Sum over all permutations of the first A; — 1 elements and over all permutations of the 
last n — k — 1 elements. 

- Sum over all possible choices for the ordered partition % = 7r 2 ). 
This is equivalent to the replacement 

EE- E E E > 

k=2 o-eS„-2 7reOP(n-2,2) creS fc _i cr'gS„_ fe _ 1 

where by OP(n — 2, 2) we denote the set of all ordered partitions of {2, 3, . . . , n — 1}. 
Furthermore, for a fixed value of k, the color factor can be written 

\ /ttia„ \ /aix \ / xa n ' v " / 

where a sum over x = 1, . . . , 8 is understood. 

Finally, the propagator clearly only depends on the choice of the ordered partition n = (jn, %%) 
and not on the order of the elements in tvi and 7r 2 . If i\\ = {n 2 , n 3 , . . . , 7r fc }, we define 

Pi,-n =Pi+ + P** + ■ ■ ■ + P*k, 

P --P { } 

1 7C2,n 1 l,7Tl ' 

At this point it is possible to identify subamplitudes in the expression for A m namely 

E 2 • • • r%fe W A k+ i ( w, . . . , a nk ,-PrX) 

crGSk-l 

= A*rl(l,1Cl-P&), (4-10) 

\ / CL y-^ \ ' 

= A_ fc+ i(-P-^,7r a ,n), (4.11) 



where x is the color of the intermediate gluon. 

Collecting all the pieces, the color-dressed BCF recursive relations read 

A n (l, 2,...,n)= E ^ ~ P i£?) An-k+i (A^^a, 

?rGOP(n-2,2) 1 ' 7ri 



(4.12) 
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Fig. 4: Diagrammatic representation of the color-dressed BCF recursive relations. 



We emphasize that although the proof of these new recursive relations relies on the adjoint 
color basis, the final result, Eq. (|4.12jl . is independent of the choice of the basis. Furthermore, 
as in the case of the Berends-Giele recursive relations, we see that the form of the color- 
dressed BCF recursive relations stays the same as in the color-ordered case, with the only 
difference that in Eq. (j4.12j) the sum goes over all partitions of {2,3, ... ,n — 1}, i.e., over 
unordered objects. This implies that the new color-dressed BCF recursive relations have the 
same properties as in the color-ordered case, namely 

1. The definition of the off-shell shifts, Eqs. (|4.2l - l4~4"|) . is independent of the color. 

2. As in the color-ordered BCF recursive relations, the pole structure of the scattering 
amplitude is manifest in Eq. (|4.12|) . 

3. Similar to the color-ordered case, the subamplitudes in Eq. ()4.12|) are not independent, 
but they are are linked via the off-shell shifts. 

The result (|4.12|) obtained for amplitudes containing only gluons can be easily extended to 
include a single quark pair. For amplitudes containing a single qq pair, the color decomposition 
reads [27] 

An(l q , 2, . . . ,n - l,n ? ) = ( T<Va • • • T^- 1 )/ 4» (l q , a 2 ,..., a n ^, n 9 ) . (4.13) 

The BCF recursive relations for color-ordered amplitudes still hold when a quark pair is 
included, where either a quark or a gluon can be chosen as the intermediate particle [17]. If 
a quark is chosen for the internal line, no sum over helicities has to be carried out, because 
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helicity is conserved all along the fermion line. The BCF recursive relations then read 

n-2 

A n (l h q , 2, . . . , n" fc ) = A k+ i (I?, 2, . . . , k, -P ff ;i) (4.14) 

fc=2 

-2— A n _ fc+1 fc + 1, . . . , n - l,fig h ) , 

q,lk 

where h is the helicity of the quark. Both the recursive relations and the color decomposition 
have the same form as in the case of a pure gluon amplitude, with the only difference that 
instead of working in the adjoint representation one now has to work in the fundamental 
representation of SU (3). So the recursive relations derived in the case of pure gluon amplitudes 
can be easily extended to include a single qq pair 

AnQ-t, 2, ■ ■ ■ , n- q h ) = A ^ *i> A n- k+ i [P h J^ n- q h ) ■ 

neOP(n-2,2) «' l7ri 

(4.15) 

The formula is exactly the same as in the pure gluon case, up to two small differences: 

- no helicity sum has to be carried out for the internal line 

- x is a color index in the fundamental representation. 

It is possible to extend this relation to include photons, by simply performing the substitution 

- tf. (4.16) 

From this it follows that a QED amplitude containing a single qq pair and n — 2 photons can 
be written in terms of color-ordered amplitudes as 

A% ED (l q ,2,...,n-l,nq)= £ A n (l q , a 2 , . . . , n q ). (4.17) 

a&S n -2 

Thus the above recursive relations ()4.14|) still hold for QED amplitudes. This particular result 
has already been pointed out by Stirling and Ozeren in Ref. [19]. However, as shown there, 
for QED processes it is more efficient to take one of the fermions and one photon as reference 
particles. In fact, as there is no photon-photon vertex, all the terms in the recursive relations 
where both fermions are in the same subamplitude vanish, simplifying the calculation. 

5 Color-dressed CSW rules 

In this section we present the color-dressing of the CSW vertex rules introduced in Refs. [6-8]. 
These rules state that it is possible to build all color-ordered amplitudes from MHV vertices, 
connecting them by scalar propagators. Each off-shell leg with momentum P then corre- 
sponds to a spinor P a a?7 a , using some arbitrary antiholomorphic reference spinor rf. However, 
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the CSW rules imply that an n-point color-ordered amplitude may have contributions from 
MHV vertices involving up to n particles. The number of different vertices is thus growing 
steadily with the number of particles, which makes it impossible to put the Berends-Giele 
recursive relations and the CSW rules on the same footing. This problem has also been ad- 
dressed in Ref. [28]. Here we introduce a new method to decompose an MHV vertex into 
three-point vertices involving an auxiliary particle and derive recursive relations for both the 
auxiliary particle and the scalar propagators. Once these relations have been obtained, the 
corresponding color dressing is performed. 

In analogy to the Berends-Giele recursive relations we define an n-point scalar off-shell 
current J^(l, . . . , n) as the sum of all MHV diagrams with n external on-shell legs, and one 
off-shell leg with helicity h. 3 This scalar off-shell current can be easily constructed employing 
the CSW rules: 



4(1, ...,n) 



pi 



n-1 



{-P{ t t p^ n ) 4,(1, . . . j) j h2 (i + 1, . . . 



n 



i=l 

n—2 n— 1 



+ 2^ ^ 4 ( ^1,"' ^iJ' ^i+lj' ^j'+l,n) • • • j : 

i=l j=i+l 

J h2 (i + 1, . . .,j)J h3 (j + l,...,n) + ... 



(5.1) 



where the dots indicate terms with higher order MHV vertices. A sum over helicities (h, h±, h,2, ■ 
with — h + hi + h<i + . . . — n — 4 is implicitly understood. According to the CSW rules, the 
vertices A n correspond to off-shell continued n-point MHV amplitudes. 

However, as mentioned above in this form the CSW relations imply a factorial growth in 
the color-dressed case, because of possibly large numbers of legs at single MHV vertices and 
the associated permutations of these legs. In order to tame this growth, we rewrite the CSW 
relations in a form similar to the Berends-Giele recursion with a tensor particle, where the 
terms in Table Q serve as basic building blocks. We first introduce auxiliary double-lines, 
carrying threefold information: 

- The total momentum P flowing in the double-line. 

- A pair (ki, k r ), describing the momenta flowing in each of the two lines separately. Notice 
that in general ki + k r ^ P. 

- A pair (a, b), describing the momenta of the negative helicity legs in the corresponding 
MHV amplitude contained in the off-shell current. If no negative helicity gluon is 
attached to the double-line, then a = b = 0, and if only one negative helicity gluon is 
attached, then b = 0. 



3 In the context of the CSW rules, it makes sense to talk about the helicity of an off-shell particle. As the 
off-shell continuation of the spinors involves an arbitrary reference spinor r] a , the scalar off-shell currents are 
not gauge invariant objects. However, the rf dependence drops out in the end [8]. 
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Right-handed Vertices 


<i i 
% J 


v K, hj ,h k( , , _ (a/3) 4 


i j 


v ag 37 k k ( l J) = S k t A r j £(hi,hj,a k b k ) 


i 3 


v aar " k k { i h % r,j) = 4 ; i,4 r j e(a i b i ,h j ,a k b k ) 


i j 


W (*^^)= <iri ) - jb ) <Wl ) e ( fl *^^ fc *) 



Tab. 1: Right-handed basic building blocks in the MHV decomposition of color-ordered 
amplitudes in the CSW approach. Details are given in the text. 

In order to build recursive relations where all n-point MHV vertices for n > 4 are decomposed 
into three-point vertices, we define n-point (off-shell) double-line currents .7^(1, . . . , n), where 
1 < u < v < n, as the sum of all diagrams with n external on-shell legs and an (off-shell) 
auxiliary double-line. This line carries the information (ki,k r ) = (Pi |U , P w +i,n) an d (a, b), a 
and b being the momenta of the negative helicity gluons attached to it. 4 Furthermore, in 
the color-ordered case, a must be of the form a = P^j with constrained to one of the 
following possibilities 

1 < i < j < u, 
u + 1 <i <j <v, (5.2) 
v + 1 <i < j < n 

and equivalently for b. Notice that according to our definition there are no one-point double- 
line currents. For later convenience we define all one-point double-line currents as zero. All 
other double-line currents can be built recursively employing the vertices given in Table ^ 
yielding 

J„t(l,...,n) = S uv V^ ab (P ltU ,P u+l!n )J hl (l,...,u)J h2 (u + l, 

v-l 

+ (1 - Suv) ^ ^AAR l '° 6 (-f , l,«! Pw+l,vPv+l,n) Juw (1> 
w=u 

where sums over repeated indices are always understood. 

4 It is easy to see that all other assignments for (ki, k r ) do not contribute in the color-ordered case. 



...,n) (5.3) 
...,v) J h (v + 1, . . . ,n), 



11 



The indices a and f3 in Table prefer to the two particles with negative helicity within one 
MHV amplitude. The e functions appearing in the vertices involving an auxiliary double-line 
keep track of the negative helicity gluons attached to them, 



1, if # negative helicity gluons < 2 
0, if # negative helicity gluons > 2 



(5.4) 



A vertex vanishes if the number of negative helicity gluons attached to a double-line is greater 
than two, since this situation corresponds to a non-MHV amplitude, that has to be decom- 
posed further by means of the CSW relations. No permutation of the incoming legs is allowed, 
since it would lead to double-counting. The respective rules in the MHV decomposition are 
obtained from the above by swaping helicities and replacing angular by square brackets. In 
analogy to the double-line current, we can write the recursive relation for the scalar off-shell 
current in terms of double-line and scalar off-shell currents 



.4(1, ...,n) 



Pi 



^ n— 1 
W 



L > n k=l 
k-1 k-1 



V^ h2 ' h (P l!k , P k+1 , n ) J hl (l,...,k)J h2 (k+l,..., n) 



(5.5) 



+ £ ^2 ^GAR h (Pl,u-, Pv+l,k, Pk+l,n) ^™(1, ■ • • , k) Jf ll (k+ 1, 



n) 



u=l v=u 



Notice that the second term vanishes for k = 1 due to the vanishing of all one-point double-line 
currents. 

The vertices given in Tabled correspond to the situation where all gluons are attached on one 
side of the double-line. We will refer to these vertices as the right-handed vertices. The right- 
handed vertices are sufficient to construct all MHV amplitudes. However, it is convenient for 
the subsequent color dressing of the CSW rules to recast Eq. (J5.5)) into a symmetric form. To 
do so, first we define left-handed vertices where all the gluons are attached to the opposite side 
of the double-line. These vertices are shown in Tabled 5 In the left-handed decomposition, 
the recursive relations Eq. ()5.3|) and Eq. (J5.5J) read 



n) 



S U v Vag 



ab (P 1:U , P u +l,n) 4i(l, ...,«) Jfca(« + 1, -. . 

+ (i - <W £ v£^> 



n) 



u 



(5.6) 



l,...,n) 



w=u+l 



n-1 



4(1, ■■■,n) = ~btY1 

1 > n k=i 



V^ h2 ' h (Pi, k , P k +i,n) 4,(1, • • • , k) J h2 (k + 1, . . . ,n) (5.7) 



n—k—l n—k—1 



^ ^ ^GAL b ' h (P^,k, Pk+l,u, P v +l,n) 4i(l, • • • , k) J!^(k + 1, . . . 



n) 



u=l v=u 



5 As Vgg and Vag are the same in both the right and left-handed decompositions, we do not list them 
again in Table [21 
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Left-handed Vertices 


I J 


V AAL 


(ijijr) = 5 hi 5 krjr e(hi,ajbj,a k b k ) 
\ l Ji) 


i 3 


* GAL 


3 ' k (%3i,jr) = ,. J,, ?»,..» eih^ajb^hk) 

{jrk){kl}{lji} 



Tab. 2: Left-handed basic building blocks in the MHV decomposition of color-ordered 
amplitudes in the CSW approach. 

Combining the right- and left-handed decompositions, it is possible to write the recursive 
relations in a symmetric form involving both right-handed and left-handed vertices 

J^(l, ...,n) = 5 uv V^ h2 ' ab (P hu , P u+1 , n ) J hl (l,...,u)J h2 (u + l,...,n) (5.8) 
1 u_1 

+ (1 " <W~ VtAR^iPhu, P W+ l, V , P V +l,n) J^'(l, ...,v)J h (v + l,...,Tl) 

w=u 

1 V t i 

+ (! - <W 5 E ^AAL ''^(^l,"' ^u+l.w, ^+l,n) 4(1, •••,«) + 1, . . . , fl), 




1 fc-1 fc-1 

+ 2 E E VS!& ,fc (^i,«, PtM-l,*, flH-l,n) JS(1, • • • , k) J hl (k + 1, . . . , n) 

u=l v=u 
- n— fc— 1 n— k— 1 

+ 2 E E V^l b ' h (Pi,k,Pk + i,u,Pv + i,n)J hl (h---,k) J2(fc + l,...,n) • 

These recursive relations, equivalent to the CSW vertex rules, can be solved simultaneously 
to construct the scalar off-shell current. The difference to the pure CSW approach without 
decomposition of the MHV vertices lies in the fact that the number of different vertices in 
Eq. (|5.8|) and Eq. (|5.9|) is fixed and does not grow with the number of particles. The approach 
thereby differs from the one presented in Ref. [28]. Furthermore, as the number of different 
vertices is fixed and as only three-point vertices are present, these new recursive relations are 
well suited to be compared to the Berends-Giele recursive relations. 

We now turn to the color dressing of the new CSW-like recursive relations. The procedure 
is very similar to the Berends-Giele case, but it contains some technical subtleties. The 
interested reader may refer to appendix A for a detailed discussion. As the new recursive 
relations only contain three-point vertices, we expect the color-dressed vertices to be of the 
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same form as in Eq. ([3.17)1 . 

V(P, Q) = 6<[6«S§ V R (P, Q) + S°6*6* V L (Q, P). 
In the color-flow basis, the scalar off-shell currents are defined in the usual way 

(1, . . . ,n) = • • • 6J r J ^n°2 • • • , <r n ). 

CT&Sn 

The color-dressed n-point double-line currents are defined by 



(5.10) 
(5.11) 



E E E 

<t'£S v - u <t"g5„_„ . 

jabf I II 

rab 



..Xir ■••^" n (5.i2) 



Tab / // / \ 

J (Tl-v)(n-u)l l <7 7r T .J 7Tm ' ^"i / 

where 717 and ir r are two proper subsets of {1,2, ... ,n} referring to the momenta flowing in 
each line separately, (A;;, k r ) = (P n , P-k t ), and 7r m is defined by 7Tj ® 7T m © 7r r = {1,2,..., n} 
(Notice that 7r m may be empty). On the right-hand side of Eq. (|5.12|) . the indices (u,v) of 
the color-ordered currents are defined by (u, v) = (#717, n — #vr r ) and 



7r m — \ n m i n m }■■■■> "mil 
7T r = {7T r ,7T r ,...,7T r ). 



(5.13) 



Finally, the symbols defined by = cr(7r|). Notice that due to the second term 

appearing in Eq. (|5.12|) . a color-dressed n-point double-line current is symmetric in (iTi,Tc r ). 
The color dressing is similar to the Berends-Giele case and the result is 



jf(l,...,n) 



PI 



7rGP(n,2) 

1 £ V^'^P.^P.^P^)^ 1 ^!,^)^^) 
7reOP(n,3) 

7reOP(n,4) 
,h\,h2,ab( n> n \<tKLi„\ tMN i 



(5.15) 



■\- ! a'b,h,ab/- D jd JD \ sra'bt ',KL(- _ \ rrMN 7 _ \ 



+ E v ^ 

7reOP(7r m ,2) 



'b',h,ab/ r> r> td \ n-a'b' ,KL( ^ „ „, \ rrMN ( _ 



+ (7Tj <-> 7T r ) 
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where in Eq. (J5.15|) OP(n m , 2) is the set of all ordered partitions of 7r m into two independent 
parts. 

Apart from a few subtleties, the procedure of the color dressing is now similar to the Berends- 
Giele case. A detailed discussion is given in appendix A. The main differences are the follow- 
ing: 

- As the color-dressed vertices, Eq. (15.10)1 . have both right and left-handed contributions, 
the symmetric form for the color-ordered recursive relations, Eq. (|5.8)) and Eq. ()5.9)) is 
employed. 

- The contributions coming from the u = v and m/o terms in Eq. ()5.8)) and Eq. (j5.9)l 
are treated separately. For example, the u = v and u ^ v contributions in Eq. ()5.9)) 
give rise to the OP(n,3) and OP(n,4) terms in Eq. (|5.14|) . respectively. 

Apart from these, the new CSW-like recursive relations retain the same form as the corre- 
sponding color-ordered relations with the difference that in the color-ordered case the sum 
goes over unordered objects. Furthermore, as in the color-ordered case, the number of differ- 
ent vertices is fixed and only three-point vertices appear in the recursive relations. Therefore 
we may compare them to the color-dressed Berends-Giele recursive relation presented in Sec- 
tion El 

6 Numerical results 

All relations for calculating multi-gluon amplitudes presented in the previous sections have 
been implemented into C++ Monte Carlo programs using the tools set AT00LS-2 . and the 
integration package PHASIC++-1 . [3]. A comparison of calculation times for helicity summed 
color-ordered amplitudes versus the results obtained in Ref. [29] has been performed. Our 
implementations yield exactly the same growth in computation time, except for the CSW 
rules, where we gain considerably due to rewriting the CSW vertex rules in terms of recursive 
relations for internal lines. Furthermore we have checked, employing the color-flow basis, that 
the color-dressed relations yield the same results as the calculations employing color-ordered 
amplitudes along with the color- flow decomposition presented in Ref. [24]. Using the adjoint 
representation, we have checked that the color-dressed BCF relations yield the same result as 
the color-ordered ones along with a decomposition of the total amplitude in the adjoint basis. 

A comparison of the computation times for the various approaches using the color-flow 
basis can be found in Table El The color-dressed Berends-Giele relations are the fastest 
method for more than five final state gluons. For less than six outgoing gluons the color- 
flow decomposition using color-ordered amplitudes calculated according to the BCF recursion 
performs better. In this case only few valid color flows exist [24] and primarily (or only) MHV 
vertices contribute. For those the computation time increases only linearly with the number 
of outgoing particles in the color-ordered BCF relations. 
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"n * i 

b mal 

State 


BG 


BCF 


CSW 


CO 


CD 


CO 


CD 


CO 


CD 




0.24 


0.28 


0.28 


0.33 


0.31 


0.26 


3g 


0.45 


0.48 


0.42 


0.51 


0.57 


0.55 


4g 


1.20 


1.04 


0.84 


1.32 


1.63 


1.75 


5g 


3.78 


2.69 


2.59 


7.26 


5.95 


5.96 


Qg 


14.2 


7.19 


11.9 


59.1 


27.8 


30.6 


7g 


58.5 


23.7 


73.6 


646 


146 


195 


8g 


276 


82.1 


597 


8690 


919 


1890 


9g 


1450 


270 


5900 


127000 


6310 


29700 


lOg 


7960 


864 


64000 




48900 





Tab. 3: Computation time (s) of the 2 — > n gluon amplitudes for 10 4 phase space 
points, sampled over helicity and color. Results are given for the color-ordered 
(CO) and the color-dressed (CD) Berends-Giele (BG), Britto-Cachazo-Feng 
(BCF) and Cachazo-Svrcek-Witten (CSW) relations. Numbers were generated 
on a 2.66 GHz Xeon™ CPU. 

It is apparent that the computation times in the color-dressed BCF and in the color-dressed 
CSW case grow very fast. In the case of the CSW relations the reason is the number of types 
of internal lines, which is larger than in the Berends-Giele and in the BCF approach. In this 
respect it is important to note that each double line may eventually carry zero, one or two 
indices of attached negative helicity gluons. Additionally, in most cases two vertices exist 
for either of these lines (cf. Table HJ, yielding a large amount of lines that finally have to be 
computed. However, the growth we encounter by employing this method is still not factorial 
but exponential. Nevertheless the factor in the exponent is still too large for the method to 
be competitive with the Berends-Giele approach. This fact is illustrated in Table 01 where 
we list the average number of nonzero internal lines counted either by value or by origination 
vertex. The former corresponds to the average number of nonzero currents in the Berends- 
Giele approach. 

Employing the color-dressed BCF relations, we encounter a factorial growth of the compu- 
tation time. We have identified three main reasons: 

- The subamplitudes are linked by the spinor shifts. 

- The natural color basis is the adjoint basis. 

- The amplitudes are decomposed down to three-point vertices. 

We address these points in order. 

In the color-dressed as well as in the color-ordered BCF relations, Eqs. ()4.1J) and (j4.12|) . 
the subamplitudes of a given decomposition are linked via the shifts Eqs. (I4.21l4.4jl . Thus 
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T~\ * 1 

.bmal 

State 


Currents 
(BG) 


Internal lines (CSW) by 


MHV vertices (BCF) by 


vertex 


value 


vertex 


value 


2g 


7.04 


3.48 


7.56 


1.98 


1.98 




19.50 


8.56 


27.45 


4.43 


4.57 


% 


44.67 


18.58 


109.0 


14.13 


18.17 


5g 


95.74 


38.63 


407.4 


63.88 


126.3 


6g 


198.8 


78.25 


1648 


297.2 


1026 


7g 


405.8 


157.8 


6773 


1395 


10330 


Sg 


850.3 


325.8 


31340 


6073 


124600 



Tab. 4: Average number of nonzero currents in the color-dressed Berends-Giele re- 
lations, average number of internal lines in the CSW approach and average 
number of nonzero MHV vertices in the color-dressed BCF relations using the 
color-flow decomposition. MHV vertices in BCF are counted either by distinct 
value or by distinct assignment of unshifted external momenta. Internal lines 
in CSW are counted either by vertex or by distinct value. 

the BCF relations need a recursive calculation of subamplitudes in the sense that the total 
amplitude is to be decomposed successively into smaller building blocks, finally yielding only 
three-point MHV vertices. In other words, we have to take Eq. (I4.12J) literally and apply 
a top-down approach of the computation, since for the evaluation of each subamplitude all 
previous spinor shifts have to be computed. Figuratively speaking, this is due to the fact that 
in the BCF recursion all subamplitudes "remember" which decomposition they originated 
from, thus inhibiting the calculation of general color-dressed subamplitudes. This fact is also 
illustrated in Table where we list the average number of distinct nonzero MHV vertices 
and the average number of distinct assignments of unshifted momenta at these vertices. The 
latter corresponds to the average number of internal lines in the CSW approach, counted by 
origination vertex. It grows much slower than the former, although faster than for example 
the average number of nonzero currents in the Berends-Giele relations. 

When applying the top-down procedure of the computation described above, it is neces- 
sary to avoid the calculation of terms yielding zero due to the color assignment of external 
and internal lines. This can be done in two steps. First, all valid color flows are identified 
employing an algorithm similar to the one used for the Berends-Giele recursion. Second, 
the subamplitudes are calculated only for the valid color structures. The calculation can be 
alleviated if the reference particles in the recursion are chosen such that together they form 
a color current having a nonvanishing contribution to the respective amplitude. Since there 
exists no decomposition assigning both particles to a common subamplitude, the correspond- 
ing color current does not contribute anymore. This procedure eliminates many terms in the 
recursion, but it is still insufficient in the case of the color-flow basis. In fact we expect some 
redundancy in the calculation of color-ordered subamplitudes due to the dual Ward identities, 
which is introduced by fixing the reference particles for all possible color flows of an amplitude 



20 



Process 


Color flow 


Adjoint 


99 ^ 2 9 


0.78 


0.86 


99 


0.83 


0.74 


99~^^9 


0.94 


0.60 


99^ % 


1.14 


0.51 


99 


1.44 


0.44 



Tab. 5: Ratio of the average number of nonzero MHV vertices in the color dressed 
and the color-ordered case in the color-flow and the adjoint representation 
decomposition. 

simultaneously, cf. Eq. (|4.12j) . To see this, consider a dual Ward identity of the form 

A(2,l,3,...,n) = - A(l,... ,1,2,1 + 1, ...,n) . (6.1) 

l=£2; l<l<n 

Assume that particles 1 and n have been fixed to be the reference particles in the recursion 
and the ordering {2, 1, 3, . . . , n} yields a valid color flow. In this case the above choice of 
reference particles is actually inconvenient to calculate the respective contribution to the total 
amplitude, since the sum on the right hand side of Eq. f|6. 1|) could be replaced by the one 
term on the left hand side. This problem does not occur in the color-ordered case, since the 
reference particles are chosen separately for each color flow. To illustrate this, in Table El we 
compare the ratio of the average number of distinct nonzero MHV vertices in the color-dressed 
and the color-ordered BCF relations for the color-flow basis and the adjoint representation 
incorporating all simplifications described above. In the adjoint representation the color- 
dressed relations yield less terms than the color-ordered ones, since the adjoint representation 
naturally avoids the problem of encountering singlet gluons, that decouple. However, much 
more effort is spent on the computation of color factors in the adjoint representation [24], 
such that it is not the method of choice. 

In the color-dressed BCF relations each amplitude is decomposed completely into three- 
point vertices. In contrast, in the color-ordered case, any MHV amplitude occuring in any 
step of the recursion can be evaluated immediately. To highlight the differences due to this 
treatment, Table El shows a comparison of the average number of distinct nonzero MHV ver- 
tices that have to be evaluated in the color-dressed and in the color-ordered case. We also 
give the same number for the color-ordered case, when each amplitude is decomposed into 
three-point vertices as well. 
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Process 


CO 


CD 


general MHV 


3-point MHV 


99 -> 2g 


1.28 


2.55 


1.98 


99 3^ 


1.84 


5.51 


4.57 


99^ ±9 


7.41 


19.33 


18.17 


99-* % 


48.78 


110.7 


126.3 


00 -> 60 


318.3 


714.7 


1026 


00^ 70 


2329 


5269 


10330 


00^ 80 


20650 


46890 


124600 



Tab. 6: Average number of nonzero MHV vertices in the color-flow decomposition for 
the color-ordered (CO) and the color-dressed (CD) BCF relations. 

7 Conclusions 

We have presented a new approach to the calculation of multi-parton amplitudes which ex- 
tends the recursive relations for the color-ordered amplitudes to relations for the full colored 
amplitudes. We have argued that in general these new color-dressed relations should be more 
suitable for a numerical implementation since they naturally avoid the factorial growth im- 
plicit in taking the sum over the permutations of the possible color flows in an amplitude. 
The taming of the factorial growth to an exponential one is easily proved in the color-dressed 
formulation of the Berends-Giele recursive relations which we find to be the same as the 
Schwinger-Dyson approach introduced in Ref. [20] and equivalent to the ALPHA algorithm 
of Ref. [30]. Using a similar approach but exploiting the adjoint color basis decomposition, 
Eq. (|2.3j) . we have proved a new formulation of the BCF relations, Eq. (j4.12jl . which involves 
the full amplitudes, including color, and retains the same formal simplicity of the original 
formulation. Finally, we have considered the CSW relations. In this case we had first to 
recast them in a form similar to the Berends-Giele relation through the introduction of a new 
type of three-point vertices and effective particles, Eqs. (j5.8l5.9J) . It is interesting to note that 
while for the Berends-Giele relations the color dressing is straightforward due to the close 
correspondence to the Feynman diagram approach, this is far less trivial for the BCF and 
CSW relations, for which there is no direct relation to the standard quantum field theory 
perturbative approach. 

To test the numerical efficiency of the different formulations we have also implemented the 
corresponding algorithms and computed squared amplitudes for 2 — > n gluon scattering, by 
performing the sum over helicities and color with a Monte Carlo method. Our results clearly 
show the numerical superiority of the recursive formulation by Berends and Giele over all 
twistor-inspired methods, both from the point of view of the growth of complexity with n 
and the simplicity of the implementation. For the color-ordered amplitude formulation we 
confirm the results of Ref. [29] except for the CSW relations, on which we improve considerably 
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by bringing them to the same level of complexity as the BCF relations. The color-dressed 
formulations of the BCF and CSW relations perform worse than the corresponding color- 
dressed Berends-Giele relations for different reasons. The BCF relations are penalized by 
their top-down structure, i.e., the fact that for each helicity and color configuration the 
decomposition in terms of amplitudes with smaller multiplicity has to be found, and from the 
fact that their natural (and minimal) color basis is the adjoint basis which is computationally 
quite heavy. The "improved" color-dressed CSW relations instead suffer from the presence 
of a large number of elementary line types and effective three-point vertices which eventually 
affect the overall growth of the algorithm. 

In conclusion, we have shown how color can be included in the color-stripped recursive rela- 
tions coming from twistor-inspired methods that do not have a straightforward relation with 
a standard perturbative Lagrangian approach. The resulting color-dressed BCF relations can 
be easily derived by employing the adjoint color decomposition, which exactly matches onto 
the BCF structure, and retain the same very elegant form of their color-ordered counterpart. 
In this respect, it is suggestive to speculate that similar colored relations might be derived 
also for one-loop amplitudes, for which an analogous color decomposition holds and may be 
deduced from an effective QCD Lagrangian, still unknown. 
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A Appendix 



In this appendix we give a detailed derivation for the color dressing of the CSW-like recursive 
relations presented in Section 03 We start with the color dressing of the recursive relations 
for the double-line current. As we want all color-dressed vertices to be of the form given in 
Eq. ([5. 10)1 . involving both right-handed and left-handed contributions, we use the symmetric 
form of the recursive relation, Eq. (|5.8jl . We will consider the terms for u = v and u ^ v 
appearing in Eq. ()5.8|) separately and we will start with the u = v term. This term corresponds 

T2l) . we obtain 



to 7r m = 0. Inserting it into the definition 



E E 

u£Su cr"£S n - u 



< a •■■^'"V^^ (A.l) 

CT r,u + l '• 1 
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The decomposition of the two color factors is similar to the Berends-Giele case, Eq. (I3.10JI . 



>/ ••• 


V ■ 


..dj 


J 

r,u-\-l 


%rT l.l 


..5 j 



Ol®N "I °i a ,i--- d K d i „ ■■■°M ' ( A - 2 ) 
<V^L 5 / °i_ „ • • • <V 5 ia,,---°K ■ ( A - 3 ) 



r,u-\-l 



As P CT7r; = P ni and P CT // = P^, we can identify two scalar subcurrents in Eq. (|A.1J) 



^fvo = E ^ (a.5) 

rr"f=<? CTr,,1+:L 

and so the u = v term reads 

vfe*"- 6 ^, po + sUM vJ&'^CP^, p^J^WJU^M. (a.6) 

We define the color-dressed AG-vertex as 6 

V^(P n , P.J = ^ 2 ' a6 (^> ^) + 4W Vjg fcl,a6 (P Wp) PrJ, (A.7) 

and so finally the u = v term reads 

V h /^ a \P^,P^ h K ; L (n l )J h A f(n r ). (A.8) 

We now turn to the u 7^ v term in Eq. (|5.8|) . This term has four contributions, corresponding 
to the right-handed and left-handed decompositions for each of the two terms appearing in 
Eq. (|5.12j) . We will show the color dressing of the right-handed decomposition of the first 
term in Eq. (j5.12j) explicitly. The three remaining contributions can be obtained in a similar 
way. The right-handed contribution to the first term in Eq. (J5.12j) gives 

E E E E^ M6 (^^^ +1 ^> P <) ( A - 9 ) 

<tGS u a'GS v - u a"£S n - v w=u 

J X^ a l>u sp"m^a X"' <T "" ja'b' 

1 . 1 

' m,u-\-l r,v-\-l 

The color factor is decomposed in the usual way 



• • • SI . . . b: ...5 j J uw {a ni ,cr n JJh{a nr ). 



XJ X^l,™ x^'m.v X^ a r," xJ xK ;M rL X^l,™ x^'m.v rJV Sp a r,r 

i(jr 1 ... i 1 =h°N<>I / i * * * it ■■■°K \„ ---V 

' m.u+1 r,v-\-l ' m.u+1 r,V+l 

(A.10) 

As P CT7r = P n and P CT // = P Wr , we can immediately identify the color-dressed scalar subcurrent 



j?"w = E € ak), (a.ii: 

^ r,u + l 



3 As for the Berends-Giele case, the color indices of the vertices are suppressed. 
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and we are left with 

E E E^f<i M6 (^p^ ro+1! ^,^) (A.i2) 

aeS u a'£S v - u w=u 

We will now consider the contributions coming from w = u and w > u separately. For w = u 
we find, with P CT / CT / = P CT / = P^ . 

6 J N Oi V AAR {P n ,P Wm ,P 7rr )d i .& i ' ...6 K J uu (a n ,a w JJ h (vr r ). 

(A.13) 

Due to the Kronecker-deltas appearing in the definition of Vaar (See Table HJ), one has 

Vit£> a \P« v P« m ,P« r ) 4U^u){< m ^n) = 0, (A.14) 
and so Eq. (|A.13|) can be written 

E E SlS^V^iP^P^P^J^M (A.15) 

cr£S u <T'(zS v - u 

( %L 5? a l,u x*°'m,v ja'b' I I \ , ri sP a 'm,v S^ a l,n ja'b' t I \\ 

[ d iv,i---\, ••• 6 K K.^ m j + \, * * iaj i ■'(e-uKti-ujC^^Tj), 

where the second term vanishes according to Eq. (jA.14j) . This second term corresponds to 
the second term in Eq. (|5.12jl 

jtM^ra) = E E -^"tfKnO (A-16) 

„ , ^ ' m,u+l 

+ \, ■■■°K J {v-u){v-u){°ir m i°n))- 

So for w = u, we get 

slsSsfv^iP^P^P^j^^in, ^ m )J h MN M- (A.17) 

For w > u, we can rewrite the sums over u> and a' in a similar way as in the Berends-Giele 
case, Eq. (j3.HI) . 

E E - E EE- ("8) 

UJ — It 

where OP(n m , 2) is the set of all ordered partitions of n m into 2 independent parts. This 
rearrangement gives, with P a n „n = P a n = P n „, 

E E E E «^f^' afe (^,^ 2 ,^) (A.i9) 

7reOP(7r m ,2) ueS u cr'€S w - u a"€.S v - w 

■■■°K JwnKiV^)^ W, 
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where n = (tti, 7r 2 ) is an ordered partition of 7r m and w = #7Ti- Due to the Kronecker-deltas 
in the definition of Vaar-, we have in a similar manner as for Eq. |04j) 



(A.20) 



Using this relation, we can identify a color-dressed double-line current in a similar way as for 
the w = u contribution. The w > u contribution reads 

E ^fv^iP^P^P^) j:% KL (n, ^W S M- (A.21) 

7reOP(7r m ,2) 

Putting together the w = u and w > u terms, we find the contribution from the right-handed 
decomposition of the first factor in Eq. (J5.12|) 



xJ xK s;M-ira'b',h,abf r> p r> \ a-db'KLl^ „, \n-MN(„\ 
d L d N d I V AAR [ P nv P 7r m ,P^) 3. ' (^h^ m )J h (7T r ) 



+ 



7rGOP(7r m ,2) 



sJ sK zM-ira'b' ,h,ab/ -n , \rra , b',KLi^ 



(A.22) 



a'b',KL/ 

7T(7T2 



rMN / 



Similar terms are obtained for the remaining three contributions. Adding up all the contri- 
butions, we find for the u 7^ v term 



a'b'.KL/ 



rMN, 



(A.23) 



■ V a A^\P^P^,P n ) j£* L (* m ,*r) J h MN (n) 

V a Z' h ' a \P^P^P,„) J a Zf\*U^) Jffrr 



E 

7rGOP(7r m ,2) 



where the color-dressed AA- vertex is defined as 

ici'b' ,h,ab / t, j~> n \ sJxKx^l j ra'b' ,h,ab / 



■>;aO,ft,oO/p p p N XJ xK cAi ytt 6 ,ft,a6/ p p p \ 1 rj rM xK T/ft,a e ,a6/ p p p \ 

"AA I-'ttd r -K 2 i r TT&) — °L°N°I v AAR > r nz) t °N°L °I v AAL \ r nz > r mi r m ) ■ 

(A.24) 
(A.25) 



Finally the recursive relations for the double-line current read 



h\,h2,ab/ 



rKLl 



rMN, 

h 2 1 



7T r 



+(1-6* 



a'b',KL/ 



rMN , 



E^ y a'b',h,ab, tj p p \ rra'b'KLf- ^ \ n-MN ( _ 
VAA i F ni , P* 2 , P Vr ) J ni n 2 ( n l > ^1 > K2J 3h (Tr 

7reOP(7r m ,2) 



+ (7T/ <H> 7T r ) 



which is Eq. f|5. 15j) stated in Section |5] 
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We now turn to the color dressing of the recursive relations for the scalar off-shell currents, 
and we use again the symmetric form of the recursive relations, Eq. (jfi.9j) . The color dressing 
of the pure gluon part is identical to the color dressing of the three-gluon vertex part for the 
Berends-Giele recursive relations, and it evaluates to 

E ^ h (P^P^ L ^)Jkf^). (A.26) 

vreP(n,2) 

The term in Eq. (|5.9jl involving a double-line current has contributions from both the right- 
handed and the left-handed decompositions. We will only show the color dressing for the 
right-handed contribution here. The left-handed contribution can be obtained in a similar 
way. The right-handed contribution reads 7 

71-1 fc-i fc-i 

EEEE^'V" ^ an ^GAR' h (P^ 1 ,a u ,Pa v+1 ,cr k ,Pa k+1 ,cr n ) Juv( a li ■ ■ ■ ) a k) Jfni&k+l, • • • , CT n )- 
crGSn k=2 u=l v=u 

(A.27) 

We consider the terms corresponding to u — v and u ^ v separately. The u = v term in 
Eq. ifPTjl reads 

71-1 fe-1 

E E E ' ' ' ^f n ^GAR' h { P n,(T u , P<r u+1 ,<T k , P<r k+1 ,a n ) J^S^U • " • > a k) ^/uO^fe+l; ■ ■ ■ , °n)- 
aeSn k=2 U=X 

(A.28) 

The color factor is decomposed in the usual way 

5^ . . . 5? - = • • • 4 fc < +1 ■ - • *M n - (A.29) 

The sums appearing in this expression can be rearranged as 

n-l fe-1 

EEE- E E E E • (a*» 

o-e5„ fc=2 «=i 7reOP(n,3) cre5 u ff'es fc _„ a"eS n _ k 

where on the right-hand side u = and k = n — #113 and OP(n, 3) is the set of all ordered 
partitions of {1, 2, . . . , n} into three independent parts. Since 

V£j_x' h (P ni , P n2 , Pn z ) J^k- U )(k- U ){^2, TTl) = 0, (A. 31) 

due to the Kronecker-deltas in the definition of Vgar, the result of the color dressing of 
Eq. (|A~28j) is 

E WC A ^,4,^)C(^^)jfN. (A.32) 

7reOP(n,3) 

Adding the left-handed contribution, we find 

E ^^(P^P^P^)^H^2)J^M, (A.33) 

7rGOP(n,3) 

7 Recall that this term has no contribution from k = 1. 
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where the color-dressed GA-vertex is defined by 

yab,h u h/p p p \ _ cJcK rM V ab,h u ht p p p \ i S.J S.M s.K v h lt ab,h, p p p n 
V GA \ r -K\i r it%i r K2,) — L°N°I 'GAR \*iri i r ^2 ) r itz) ^ °N°L °I v GAL \ r n 3 1 r tti j r -K%) ■ 

(A.34) 

The color dressing of the u 7^ v contribution in Eq. (JA.27J) is similar to the u = v case, except 
that we have to introduce the set OP(n, 4) of all ordered partitions of {1, ...,n} into four 
independent parts, and the rearrangement of the sums is now written 

n-l fe-1 fc-1 

EEZE- E E E E E • (aj») 

creS„ fc=l u=l v=u+l w£OP(n,4) a£S u a'€S v ^. u a"£S k - v a'"eS n - k 

where on the right hand side u = #tti, v = j^T\\ + #7r 2 and k = n — #7T 4 . The result of the 
color dressing of the u ^ v contribution then reads 

E V^IP^P^Fj^Vi.^.^Uf^). (A.36) 

7rGOP(n,4) 

Finally the color-dressed recursive relations for the scalar off-shell current read 



p2 



E ^^,^ 2 )^f (A.37) 

vreP(n,2) 

1 ^ vS^l^.^.pg^Vi.^Cfe) 



2 

7rSOP(n,3) 



+ 5 E ^'^p.^p^^.J^f^K,^,^)^^) 

7reOP(n,4) 

which is Eq. (j5.14|) stated in Sectional 
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